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Abstract 

We use the discrete kinetic theory with the free-orientation parameter being fixed 
to derive the macroscopic velocity field for many particles flowing through a microdomain. 
Our results resemble qualitatively other hydrodynamical solutions. The V-shaped velocity 
field changes as the dominant physical parameter (Knudsen number) varies. We also briefly 
discuss the possible mechanism due to the entropy production along the boundaries. 
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1 Introduction 

Discrete kinetic theory [1-2] with the thermostat assumption or diffuse scattering boundary con- 
dition [3] have been adopted to ad hoc model the many-particle scattering situation along the 
gas-surface interacting interface in a plane channel [1] . Specific orientations relevant to different 
rarefaction parameters were identified therein [1]. Motivated by the recent interests in the quan- 
tum Boltzmann approach [4] and the relevant studies [5-8], we continue our previous studies [1] 
by examining the related velocity and vorticity fields corresponding to those specific 9 and Kn 
we obtained and checking the special case 8 = 7r/4. 

Many interesting problems have been successfully solved [9-17] by using the discrete kinetic the- 
ory. Carleman (1957) developed 2-velocity models which are defined by abstract properties in 
order to produce only mathematical purposes. That model was not constructed on the basis of 
detailed analysis of the collision mechanics [9]. Broadwell (1964) devised a 6- velocity model to 
handle the simple shock wave problem [10]. At first mathematicians have been interested in the 
discrete Boltzmann equation with the hope to provide, thanks to the relatively simpler struc- 
ture of the equation as compared with that of the full Boltzmann equation, stronger results than 
those available for the Boltzmann equation or mathematical results suitable to be technically 
generalized to the full Boltzmann equation in the analysis of initial and initial-boundary value 
problems. However, the analysis over recent years has shown that this is not the case [5,7-8]. 
These have been reviewed considering mainly the mathematical aspects of the initial and initial- 
boundary value problems in order to provide a "more detailed" analysis in a "more specialized" 
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field. In fact the consistency of the mathematical results depends on the structure of the math- 
ematical problems : in some cases it is possible to obtain for the discrete Boltzmann equation 
"stronger" results than the corresponding ones for the Boltzmann equation, and in other case 
"weaker" results. Kawashima has proved the global existence of solutions to the initial-boundary 
value problems (I.-B.V.P.) for the discrete Boltzmann Equation in the ID-region < x < oo or 
< x < d (cf [11]). 

Cornille obtained some transient or stationary family of solutions for certain (fixed-orientation) 
discrete velocity models without considering the boundary conditions [12]. 
In this short paper, we plan to reconstruct the macroscopic velocity field of many dilute particles 
by the verified 4- velocity model [18] (the free orientation is fixed to be 8 = vr/4) considering a 
simple test problem : many molecules or particles flowing along the bounded-plane channel and 
finally reaching a steady state. The verification of our approaches has been done in [1,6] (cf. 
Chu), the argues about the differences between different discrete velocity models included. For 
= 7r/4 case, using a completely different solving procedure, we obtained velocity fields which 
have a V-shaped or chevron-structure. 

This short note is organized as follows. We introduce the general orientation-free [5,6] 4-velocity 
model in Section 2, and simplify it to a system of four equations for associated unknown func- 
tions. The general boundary conditions will be briefly introduced, too. Then, we define some 
macroscopic variables (like u, v) to suit our interest which is to find a class of steady (and paral- 
lel) non-boundary-driven solutions or flows for particles flowing along a microslab with bounded 
(flat-plane) walls. The orientation will be fixed as 7r/4 here when we solve the time- independent 
system of equations with relevant boundary conditions for the test case. As reported in [18], 
there will be no dispersion or absorption when we implement the model with this orientation so 
that we can resolve sharp velocity profiles. These kinds of solutions, u, which collect the main 
results of the present paper, are given in explicit form, and are functions of ID coordinate : y 
or Y and are also dependent on certain integration constants or parameters due to the purely 
diffuse reflection boundary conditions. Finally, we analyze the solutions (V-shaped fields) and 
make some physical comments or discussions in comparison with the other flow-pattern selection 
due to the relative orientation effect upon binary encounter of many particles or unusual entropy 
production along the confined boundaries. 

2 Formulations 

Considering a simple fluid of particles with mass m and cross-sectional area a, the first step 
of the modelling procedures consists in discretizing the velocity directions in a finite number of 
unit vectors ik, k = 1, • • • ,p. One or more moduli are then associated to each direction. The 
ratio of the moduli has, however, to be properly chosen, so that collisions between particles 
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with different velocity moduli are possible. For one velocity moduli case, Uj= cifc, k = 1, ■ ■ ■ ,p; 
c = c(x, t) in general. Normally c is determined by the equilibrium distribution. 
The particles (hard-sphere) move in the whole space and collide by simple elastic collisions 
locally in space. The mathematical model is an evolutional equation for the number densities 
iVj(x,i) linked to the finite set of velocities Uj. We write a balance equation for the number 
density of particles " i" in the form 

[-| + m ■ V]N t = Gi-Li 

where Lj and Gi are the loss and the gain of the particles "i" due to collisions. In case of binary 
collisions an exact balance may be obtained, and is expressed with the transitional probabilities 
and the number densities. This model has the structure of a system of semi- linear partial 
differential equations of hyperbolic type. Above equation could also be written as 

-Ni + u i -VAT i = E E E S(i,J r ,I r )A J / r N Ir , 

r=2I r <=E r J r £E r 

where i = 1, ■ ■ ■ ,p; here, by definition, an r-collision (r > 2) involves r particles. I r =(i\, ■ ■ ■ , i r ), 
and J r =(ji, • • • ,j r ) are two elements of E r , which is the set of r-not arranged numbers (con- 
sidering the combinations instead of the order they appear) taken in the set {1, ■ ■ ■ ,p}. 
A " transitional" probability denoted by A 3 f is associated to each r-collision I r — ► J r . In the case 
of binary collisions, this term (also is called as the transition rates) is referred to the collisions 
(uj,Uj) <-> (iifc,uj), i,j,k,l = 1, and the number of paired-outputs corresponding to a 

given paired-input is denoted by q. Ni r denotes the product N^N^ • • • N-i r . 
5(i, J r , I r )= 5(i,J r )— 8(i,I r ) is the algebraic number of particles "i" created through the colli- 
sion I r — > J r . 5(i,I r ) is (positive or zero) the number of indices i present in the r-set. If only 
nonlinear binary collisions are considered and considering the evolution of iVj, we have 
r) N- v 

3=1 (k,l) 

where (k, I) are admissible sets of collisions. We may then define the right-hand-side of above 
equation as 

3,k,l 

with i £ A ={1, • • • ,p}, and the summation is taken over all j, k, I £ A, where A l l, are nonneg- 
ative constants satisfying A^ = A 1 ^ = A\ 3 k : indistinguishability of the particles in collision, 
A l ^(ui + Uj — Uk — ui) = : conservation of momentum in collision, A 1 ^ = A™ : microre- 
versibility condition. The conditions defined for the discrete velocity above requires that elas- 
tic, binary collisions, such that momentum and energy are preserved Uj + Uj = Uj, + u/, 
|uj| 2 + |uj| 2 = |ufc| 2 + |u^| 2 , are possible for 1 < i,j, k,l < p. 
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The collision operator is now simply obtained by joining A^j to the corresponding transition 
probability densities a^j through A*j =S\ui — uj\ afj, where, 

a*j>0, ]T a $ = 1 > Vi,j = l, 

fc,/=i 

with S 1 being the effective collisional cross-section. If all q (p = 2q) outputs are assumed to be 
equally probable, then a^J=l/g for all k and /, otherwise a^j= 0. The term S\\ii — Uj\dt is the 
volume spanned by the particle with Uj in the relative motion w.r.t. the particle with Uj in the 
time interval dt. Therefore, S\ui —\ij\Nj is the number of j-particles involved by the collision 
in unit time. Collisions which satisfy the conservation and reversibility conditions which have 
been stated above are defined as admissible collisions. 
The discrete kinetic equations then [1,5,7,18] assume the following form 

^1 + c[cos[e + (<-!)* n/q)^ + sin(# + (i - 1) * 7r/g)S = ^£ (N,N, +q - 

7=1 ' 

2cS 9-1 

iViiV i+g ) or = ^(Ni+iNi+i+g-NiNt+q); % = l,---,2q, (1) 

9 1=1 

where 9 is the free orientation starting from the positive x— axis to the u\ direction [1,18], 
Ni = Ni + 2 q are unknown functions, and c is a reference velocity modulus. 

According to [13], for the 2q- velocity model that is q > 3, there are more collision invariants 
than the physical ones or conservation laws which are corresponding to the number of macro- 
scopic variables (in 2D. there are only 4, i.e., one mass, two momenta, one energy). That's to 
say, there are unphysical or spurious invariants or macroscopic variables for q > 3 models (which 
could be, however, well handled by adding multiple collisions [13]). Thus, we plan to use only 
the orientation-free 4-velocity model for our test-case problem below. 




Fig. 1 Reference frame for the 4-velocity model with 9 = ?t/4 here. 
2.1 Boundary Conditions 

We use purely diffuse reflection boundary condition [1,3,15-16] here, which means properties of 
the reflected particles are independent of their properties before the impact. In other words, the 
re-emitted stream has completely lost its memory of the incoming stream, except for the conser- 
vation of the number of particles. Moreover, we impose the following conditions : the particles 
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are in Maxwellian equilibrium with the wall ("the wall locally behaves as a thermostat", i.e., the 
particles reflect after they have been in thermodynamic equilibrium with the wall-temperature) 
satisfies Ni(r,t)=ji(r,t)N w i(r,t), where 7$ expresses the accomodation of the particles to the 
wall quantities, and N w i is the discrete Maxwellian densities for the 'i'-direction set of particles. 
That is, we have 

|u_j • n\N wj = Bij\ui ■ n\N wi , j G R, B { j > 0, Y B n = ^ ( 2 ) 

iei j&R 

with / = {i, (uj — u w ) n < 0} related to the impinging set of particles, R = {j, (uj — u w ) n 
> 0} related to the emerging set of particles, n is the outer normal, u„ is the wall velocity. 

2.2 Flows in a Plane Channel 

We firstly define the related macroscopic variables n = N1 + N2 + N3 + N4, nil =c{aN\ — (3N2 — 
aN^ + PN4), nV = c{(3N\ + 01N2 -PN3 — aN^), (the latter two are the momentum flux along 
x- and y-directions) with p = n m, m is the mass of the molecule, p is the macroscopic density 
of the gas. Then, set rij = Ni/n, i = 1,2,3,4; and then use non-dimensional u = U/c, v = V/c, 
Y = y/d, where c may be related to the external forcing [1,19]. d is the full channel width. 
y = is along the center-line. 

The geometry of a 2D problem we shall consider is a kind of microchannels with bounded flat- 
plane walls which are separated apart by a width d. Particles (driven by an external constant 
forcing initially) flowing along this channel will finally reach a fully developed state (steady state 
and du/dx = 0, v = 0). 

We derive the solutions with a(= cos 9) = /?(= sin#) = \/2/2 case here. The algorithm is differ- 
ent from those previously reported, we must solve the independent number density respectively 
then combine them into macroscopic ones since the original macroscopic equation is singular 
(cf. equations in [6] by Chu). Meanwhile, from the preliminary results reported in Ref. [1,18], it 
seems, for the case of 6 = vr/4, 4- velocity model will give completely different dispersion relations 
for the thermodynamic checking of the perturbed Maxwellian equilibrium state. There will be 
no dispersion or absorption for this particular case. 

The governing equations (1), for the assumptions prescribed above, now become 
dni dn 2 dn 3 dn^ V2 ^2 

here, 713 = a — ni, n 2 = 1/2 — n±, = 1/2 — a + rii; Kn= l/(dSn) is the Knudsen number. The 
diffuse reflection boundary conditions become : 

N w2 N 1 = N wl N 2 , PN X + aN 2 - 0N 3 - aN 4 = 0, (4) 

it means (i) the Maxwellian equilibrium at the walls dominates, (ii) no penetration occurs across 
the wall. The discrete Maxwellian densities N w i at the wall, as derived before (please see the 
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detailed references in Refs. [14-16]), are 

N wi = (n/4){l + (2/c 2 )u w • u, + (-1)*[K • u 2 ) 2 - (u w • Ul ) 2 ](l/c 4 )}. (5) 

Here, boundary conditions = (the walls are static and fixed) and by assuming the 

symmetry, 

ni = B 31 n 3 + B 4 in 4 , at Y = -1/2, n 3 = B 43 ni + B 23 n 2 , at Y = 1/2, (6) 
with the discrete Maxwellians N w i\± = 1/4. Integration of Eq. (3) gives 

4Kn 

Now, set A = l/(4\/2 Kn), so we get from above equations to solve for a, b : 

[2A{1 + B 31 - B 41 ) - B 31 + B 41 ]a + (1 + B 31 - B 41 )b = -^ + {l + B 31 - B 41 )A, (7) 

[1 + 2A(l + B 13 -B 23 )]a + (B 23 -B 13 -1)6 = ^ + A(l + B 13 -B 23 ), or (8) 
Ca + Db = G, E a + F b = H. (9) 

After manipulations, we have 

FG-DH , CH-EG . , 

a= CF^DE> b= CF-DE ] (10) 

where CF - DE = AA{B 23 - B 13 - 1) + (5 23 - B 13 ){B 41 - B 31 ) - 1, and FG - DH = 

(B 23 -B 13 -l)[A{l + B 31 -B 41 )+B 41 /2] + {l+B 31 - B 41 )[A{B 23 -B 13 -l)-B 23 /2], CH-EG = 

A[(B 23 - 1)(1 - B 4 i) - £ 3 i(2 + £ 13 - 2£ 31 )]+[(£ 23 - 1)^41 - B 23 B 3l }/2. 

Since n£7/c = y/2/2(N\ — N 2 — N 3 + A4), so we have a family of (particular) flow field in terms 
of the macroscopic velocity 

u=V2(2n 1 -a) = ^^Y + V2(2b-a). (11) 

3 Results and Discussions 

This class of solution u obtained by fixing the orientation to be 7r/4 is in general different from 

those reported in Ref. [6] by Chu. Note that, for one extreme case of boundary conditions as 

mentioned in Eq. (2) : B 34 = B 4 \, and B 43 = B 23 ; we have 

_ 2A + (£41 + B 23 )/2 u _ A (l + B 41 -B 23 ) + B 41 /2 

a ~ 1 + 4A ' b ~ 1 + 4A ■ [ ] 

We can easily observe that, from equation (11), that u = everywhere for all Knudsen numbers 

(Kn). There is no macroscopic flow [20] for many particles once the boundary conditions are 

selected above. 

Otherwise, the velocity field (from equation (11)) as shown in figure 2 is qualitatively similar to 
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the V-shaped or chevron-like structure or pattern [21] reported before in other physical systems. 
The velocity field is tuned mainly by the Kn and weakly by a and b with the latter due to 
the boundary conditions. We note that a might depend on the physical properties of fluids 
and the geometry of the solid-wall as it comes from the gas-solid interaction or reflection. The 
flow-pattern selection mechanism is yet open to the best knowledge of the authors but might be 
partially linked to that reported in Ref. [1,6] (cf. Chu) since there will be an essential singularity 
when integrating equation (1) for 6 = ir/A case. In short, as Kn increases, the chevron front 
becomes more flat. 

The macroscopic vorticity to (or the mean shear) could be obtained by noting 

du I -2a d[V2(2b-a)] 
UJ -dY-~K^ + dY ' (13) 

with 

= 2A(1 + B 31 - B 41 )(B 23 - gig - 1) + [2B 41 B 23 - B 41 (B 13 + 1) - B 23 (l + B 31 )]/2 
° 4A(B 23 - B 13 - 1) + (#23 - B 13 ) (B 41 - 531) - 1 

where the last term of equation (13) is generally zero. Once the Knudsen number (Kn; a kind 
of rarefaction measure for many-particles interactions or collisions) is fixed, the vorticity is a 
constant with the related B^ subjected to the constraint in equation (2). In fact, B^ should 
depend on the detailed interactions of the gas-solid interface, like a kind of (known) molecules 
colliding with specific walls made of (already) specified material. It is bounded above but diffi- 
cult to be fixed even for specific model and boundary value problem [11]. 

Our results for the vorticity field, at least, qualitatively matches with the hydrodynamic two- 
dimensional solution [20] when the weakly compressible (incompressible) particles flow along a 
static flat-plane channel and finally reach a fully developed state even though the particles are 
initially driven by a constant pressure-gradient or unit forcing. Interestingly, similar sharp flow 
fields of solitary wave profiles (the highest one, cf. Figs. 9 and 10 by Wu et al. in [22]) and 
constant-V vortex was reported recently [22] in other physical system dealing with confined flow 
transports. 

To further interpret the mechanism, we propose that the complicated rate of entropy produc- 
tion along the boundaries (cf. [23-24]) might favor the smearing of viscous diffusion (toward the 
away-from-wall regions) so that the sharp and strange pattern could form and then there is no 
significant smoothing of the profiles along the cross-section. 

Note that, the approach here : firstly tracing or obtaining (solving the corresponding equation 
in (3)) each individual discrete number density (n^) then by summing up the corresponding 
projection to obtain u, is different from that in [6] (by Chu) : directly construct the macroscopic 
solutions from the relevant governing equation for macroscopic variables (u). The boundary 
treatment which is relevant to the entropy production there is thus entirely different. The cor- 
responding non-equilibrium states (due to different rates of entropy production and their decay) 
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approaching to the final equilibrium states which are used as our boundary conditions might 
then be different. In fact, as we noticed, the argument raised in [24] could be applicable to 
present approach as evidenced in the boundary operator as expressed in equation (2) (could 
be represented as similar divergence form). Otherwise, if our interpretations don't work, there 
might exist other unknown mechanism which need our further works. 
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Fig. 2 (a),(b) Rarefaction effects (Kn) on the velocity field u or the 
V-shaped or chevron-like structure. Kn= l/(dSn) is the Knudsen number. 
S is the effective collision cross-section, n is the number density of particles. 



